APS/123-QED 



First-Order Quantum Phase Transition in Bose Gases 

Nguyen Thanh Phuc/ Yuki Kawaguchi,^ and Masahito Ueda^ 

^Department of Physics, University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 
^Department of Applied Physics and Quantum-Phase Electronics Center, 
University of Tokyo, 2-11-6 Yayoi, Bunkyo-ku, Tokyo 113-0032, Japan 
(Dated: January 17, 2013) 

The BogoUubov theory has proven to be an accurate and versatile tool in the study of weakly- 
interacting dilute Bose gases at low temperatures, yet there is one exception where it goes qualita- 
tively wrong, i.e., near the first-order quantum phase transitions. By examining a phase transition 
in spinor Bose-Einstein condensates (BECs), we find that the energy spectrum given by the Bogoli- 
ubov theory is inconsistent with the fact that the phase transition is first order. We resolve this 
problem by calculating the spectrum based on the spinor version of the Beliaev theory. We also 
discuss the ground-state phase diagram of spin- 2 BECs which is modified by quantum fluctuations 
and the possibility of macroscopic quantum tunneling near the cyclic-nematic phase boundary. 
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Quantum phase transitions have been an active field of 
research [l|, Q , and played a central role in the study of 
electrical and magnetic properties of solid-state materi- 
als such as magnetic insulators, heavy fermion materials, 
novel semiconductor structures, and high- temperature 
superconductors. In ultracold atoms, the quantum phase 
transition from superfluids to Mott insulators has been 
extensively investigated both theoretically and experi- 
mentally jsl, 3 . Although the focus of these studies has 
been mainly on the second-order or continuous quantum 
phase transitions, the first-order quantum phase tran- 
sitions in fermionic systems such as itinerant-electron 
magnets [H and superfluid helium 3 [6[ have recently at- 
tracted much attention in connection with an extended 
non- Fermi liquid phase and superconductivity 0, HI . 

In this Letter, we point out that the Bogoliubov the- 
ory jsl , which has been widely used to investigate weakly- 
interacting dilute Bose gases [10], is inadequate to de- 
scribe spinor Bose-Einstein condensates (BECs). In 
spinor BECs, there exist several ground-state phases with 
distinct symmetries [11]. The phase transitions between 
them, therefore, should be first order. However, the re- 
sult of a Bogoliubov analysis implies that there is no 
metastable state around the phase boundary, suggesting 
that the phase transition is second order. We resolve this 
inconsistency by invoking the spinor version of the Be- 
liaev theory 12|, ll3| . We will elucidate the above state- 
ment by taking spin-2 Bose gases as an example. Af- 
ter showing that quantum fluctuations modify the mean- 
field phase diagram of spin-2 BECs, we examine in detail 
the transition between the uniaxial nematic and cyclic 
phases and prove that the Beliaev energy spectrum ac- 
tually supports metastable states. We argue that this 
should be a generic feature of first-order quantum phase 
transitions in bosonic systems since the energy landscape 
of a first-order quantum phase transition requires higher- 
order terms that are beyond the (j)^ model and Bogoli- 
ubov theory. Finally, the possibility of quantum macro- 



scopic tunneling (MQT) from a metastable state to the 
ground state in the first-order quantum phase transition 
is discussed. 

The Hamiltonian of a homogeneous spin-2 spinor BEC 
is given by H = ho -\- V with kinetic energy ho = 
fi^V^/2M)'0,(r) (M is the atomic 
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Here, :: denotes normal ordering of operators, and h = 
E,-^](r)^i(r), F^= E.jV'J(r)(f)iii-(r), and ioo = 
{l/^/b)^J{—l)~^^lJj{r)ljJ-j{r) are the number density, 
spin density, and spin-singlet pair amplitude operators, 
respectively, with (f)^^ being the spin-2 matrix vector. 
The coefficients co,ci,C2 are related to the s-wave scat- 
tering lengths ajr (J-* = 0, 2, 4) of the total spin channel 
by Co = 47rn2(4a2+3a4)/(7M), a = ATTh^{a^-a2) / [IM), 
C2 = 47rfi^(7ao - 10a2 + 3a4)/(7M), respectively. The or- 
der parameter is represented by a five-component spinor 
$ = V^fe,6,^'o,^-i,^-2)^, where no is the number 
density of condensed atoms and T denotes transpose. 
At the mean-field level, the ground-state phase diagram 
obtained by minimizing the energy functional is shown 
in Fig. [TJa) [l^. There is an energy degeneracy in 
the manifold of nematic phases with order parameters 
= y^(sin 77/^/2, 0, COST], 0, sin 7^/ v/2)^, where 77 is the 
degeneracy parameter. 

Since the order parameters of different phases cannot 
be transformed continuously to each other, the phase 
transitions between them must be first order. Neverthe- 
less, the energy spectrum given by the Bogoliubov theory 
shows either a dynamical or a Landau instability at all 
phase boundaries of spinor BECs under zero magnetic 
field. For example, one of the elementary excitations 
of the nematic phase has the energy spectrum given by 
hujp = ^eO(eO - 2c2n), where = h'^p'^/{2M) and n is 



2 



the total number density. This spectrum becomes a com- 
plex number at long wavelengths for positive C2 , indicat- 
ing a dynamical instability of the nematic phase. Such an 
instability at the phase boundary, however, characterizes 
second-order as opposed to first-order pha se transitions 
with the existence of metastable states [15|. The Bogoli- 
ubov energy spectrum is, therefore, inconsistent with the 
fact that the phase transition is first order. 
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FIG. 1: Ground-state phase diagram of spin-2 BECs ob- 
tained with (a) the mean-field approximation and (b) the Lee- 
Huang- Yang (LHY) correction. The dashed lines in (b) indi- 
cate the mean- field phase boundaries in (a). The inset in each 
phase shows the surface plot of the spinor order parameter de- 
fined as = I X^m where the hue indi- 
cates the phase of 0) according to the color gauge shown 
on the right. The LHY corrections due to quantum fluctua- 
tions lift the degeneracy in the nematic phase in (a), leading to 
the uniaxial nematic (UN) and biaxial nematic (BN) phases 
for ci > and ci < 0, respectively. Quantum fluctuations 
also shift the cyclic-UN and ferromagnetic-BN phase bound- 
aries. However, the cyclic- ferromagnetic phase boundary is 
not affected by quantum fluctuations because of the destruc- 
tive interference caused by time-reversal symmetry breaking 
(see text). 

We first show that the mean-field phase diagram is 
modified due to the Lee-Huang- Yang (LHY) correction, 
which arises from quantum fluctuations. As shown in 
Refs. [l6|, ll7|, quantum fluctuations lift the energy de- 



generacy in the manifold of nematic phases, making the 
uniaxial nematic (UN) phase (77 = 0) and biaxial nematic 
(BN) phase (if] = 7r/2) the ground states for ci > and 
ci < 0, respectively. Quantum fluctuations also shift the 
UN-cyclic and BN-ferromagnetic phase boundaries as fol- 
lows. Up to the LHY correction, the ground-state energy 
densities of the cyclic and UN phases for ci > 0, C2 < 
are given by [l8[ 
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respectively. Here, we expanded E^^ in powers of 
C2/C0 and C2/C1, which are expected to be small near 
the UN-cyclic phase boundary. In fact, since \c2\ ~ 
^3/2^1/2^5/2^^3 p^^g^ boundary [see Eq. g])], the 

last term in Eq. (|3]) is smaller than the others by a factor 
of vW^ <c 1 with a = (4a2 + 3a4)/7 = coM/(47rfi^), and 
thus can be neglected. By comparing the energies in Eqs. 
(j2j) and (|3j), we find that the phase boundary between the 
UN and cyclic phases is shifted to 
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i.e., the region of the cyclic phase is expanded. Similarly, 
the ground-state energies of the ferromagnetic and BN 
phases with LHY corrections for ci < 0, C2 — 20ci are 
given by [l8| 
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respectively. By noting that |c2 — 20ci| ~ 
^3/2^1/2 1^^|5/2/^3 ^g^^ ^Yie phase boundary [Eq. Q], 

the last term in Eq. (|6]) is smaller than the others by a 
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factor of Vna^ <C 1, and thus negligible. Consequently, 
the phase boundary between the ferromagnetic and BN 
phases is shifted from its mean-field value C2 = 20ci to 
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i.e., the region of the ferromagnetic phase is expanded. 

On the other hand, the phase boundary between the 
ferromagnetic and cyclic phases at ci = and C2 > is 
not influenced by quantum fluctuations since both of the 
ground-state energies are independent of C2, and equal to 
each other at ci = [see Eqs. (|2]) and (|5])]. That and 
are independent of C2 is due to time-reversal symme- 
try breaking in these phases. For the ferromagnetic phase 
with order parameter (j) = y^(l, 0, 0, 0, 0)^, excitations 
caused by interaction C2 such as (p = 0, = 2) + (p = 
O^rnp = —2) (p,mF = 1) + (— p,mF = —1) do not 
occur due to the absence of condensed atoms in the mag- 
netic sublevel mp = —2. Meanwhile, the cyclic phase 
has the order parameter (j) = y^(l/2, 0, a/2/2, 0, —1/2)^ 
with opposite sign for = ±2 components, which 
results from time-reversal symmetry breaking. There- 
fore, the amplitudes of excitations from (p = O^mp = 
2) + (p = 0, = —2) and from (p = 0, uif = 0) + (p = 
0,mi? = 0) to (p,mi? = 1) + (— p,mi? = —1) both via 
interaction C2 interfere destructively with each other. 

In the presence of an external magnetic field, the 
difference in the LHY correction among ground-state 
phases, whose magnitude is of the order of AE 



Green's function-based Beliaev theory is applied to cal- 
culate the second-order energy spectrum of the UN phase 
in a manner similar to our previous work on spin-1 spinor 
BECs The Dyson's equation for the Green's func- 
tions is given by 

where the interacting Green's functions, non-interacting 
Green's functions, and self-energies are denoted by G, 
G^, and E, respectively, all of which are 10 x 10 ma- 
trices with j,j',m,m' = — 2, . . . , 2 indicating the spin 
index and a, /3, 7, J = 1,2 accounting for the normal and 
anomalous components. For the UN phase, each of the 
10 X 10 matrices reduces to five 2x2 submatrices, which 
facilitates the analytic calculation. The zero-momentum 
energy of the elementary excitation that drives the tran- 
sition from the UN to cyclic phase is then calculated to 
be 
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where /u^^^ is the second-order chemical potential, and 
A, B, and C are the lowest-order coefficients in the 
frequency-power expansion of the second-order self- 
energies T,l]^'^\p), T,l^P{p), and T,l^^_^^{p), respectively. 
These coefficients are given by 
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M^/^Ci^^n^/V^r^^^^ can compete with the quadratic Zee- C Vl\fZc-l^ (^^^0 - 30v3ci j C2 



man energy g^, and the phase diagram depends on the 
relative strength of these two effects. In the limit of a 
high magnetic field A£^, the effect of quantum fluc- 

tuations can be neglected, and the ground-state phase 
diagram is obtained by mean- field theory [l9|. In the 
opposite limit of low magnetic field qB ^ A^, the quan- 
tum fluctuations dominate, the quadratic Zeeman energy 
becomes negligible, and the ground- state phase diagram 
would be Fig. [IJb). These two distinct regimes can be 
crossed over at ~ which corresponds to a mag- 
netic field of the order of 7 mG by using the parameters 
of ^'^Rb 20|,|2lj and atomic density n = 10^^ cm~^. If an 
external magnetic field is made as small as 0.1 mG (22j . 
all the above regimes can be investigated experimentally. 

The inconsistency between the energy spectrum given 
by the Bogoliubov theory and the phase diagram even 
becomes more serious with the addition of LHY correc- 
tion to the ground-state energy. Indeed, the cyclic phase 
is the ground state for c™~^ < C2 < whereas the 
Bogoliubov energy spectrum shows an instability of the 
cyclic phase for C2 < 0. Now we show that all the above 
inconsistencies are resolved if we go to the next level of 
approximation, i.e., with the Beliaev theory 12, 13|. The 
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Here, using the same argument as below Eq. ([3j), all 
terms with higher powers of C2 are smaller than those in 
Eqs. (p!Q]) - (p!2]) by at least a factor of vW^ <C 1, and thus 
neglected. By substituting Eqs. (p!Q]) - (p!2]) into Eq. (|9]), we 
find that the UN phase becomes dynamically unstable if 
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From Eqs. (j4]) and ([T3j), we see that the UN phase be- 
comes metastable for 
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The existence of this metastable state and the associated 
hysteresis phenomenon are characteristic of first-order 
phase transitions [l5|. Similarly, a metastable state of 
the cyclic phase is expected to exist in the parameter 
regime c™"^ > C2 > c^"^"'*^^^^ Consequently, the in- 
consistency of the Bogoliubov energy spectrum with the 
first-order phase transition and the ground-state phase 
diagram has been resolved by using the Behaev theory. 

We notice that such a disappearance of metastable 
states should be a generic feature of the Bogoliubov 
theory for any bosonic system. The Bogoliubov the- 
ory is included by the dynamical form of the mean- 
field theory characterized by the time-dependent Gross- 
Pitaevskii equation [lo|- At this level of approximation, 
the field operator is replaced by the order parameter, i.e., 
the condensate wave function ^, and the energy func- 
tional is described by a (j)^ model. On the other hand, in 
the case of a first-order phase transition, terms beyond 
(j)^ become relevant [l5|, and they can only be obtained 
by using models such as the Beliaev theory that take into 
account quantum fluctuations. 

To make a quantitative comparison with other first- 
order phase transitions, we now evaluate the size of the 
first-order quantum phase transition between UN and 
cyclic phases. For a first-order thermal phase transition, 
its size is defined by the ratio of the latent heat to the 
jump in the specific heat and shown to be proportional 
to the size of the temperature regime of the metastable 
state, i.e., the difference between the transition temper- 
ature and the temperature T^"^*^^^® at which the in- 
stability occurs [2^. Renormalized by T^, the dimen- 
sionless size S = |(T^ — T^unstabiej^^^cj normal to 

superconducting transition is S ~ 10~^, and that for the 
nematic to smectic-A phase transition in liquid crystals 
is 6 ^ 10~^. This indicates that they are weak first-order 
phase transitions induced by thermal fluctuations 0. 
Similarly, the dimensionless size of the first-order quan- 
tum phase transition under consideration can be evalu- 
ated as 



^UN-C _ ^UN-unstable 



UN-C 
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Using Eqs. dU and ([T3|), we obtain S 0.1. This value 
is much larger than the corresponding values of S for 
the above thermal phase transitions, indicating that the 
quantum phase transition between the UN and cyclic 
phases with distinct order parameters is a rather strong 
first-order phase transition. 



From Eq. ()T4)) . there is a parameter regime in which 
the UN phase is metastable while the cyclic phase is the 
ground state. By neglecting a small quantum depletion. 



these states are approximately given by 
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where aj^ is the creation operator of a particle with zero 
momentum and magnetic quantum number mi? = m for 
the vacuum state |vac). Because these states are not ex- 
act eigenstates of the many-body Hamiltonian ([B , they 
can undergo quantum diffusions in the spin space |25l-[28|. 
and induce a MQT. We now estimate the time scale of 
the MQT by using a two-state model, i.e., by restricting 
the Hilbert space to these two states with local energy 
minima. The time scale of the MQT is then given by 
r = h/A with A = 2(Cyclic|T>|UN). Using Eqs. ([IS]) and 
([1]), we obtain 



2N/2 



con(iV-l)' 
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where N is the total number of particles. The exponen- 
tially large factor of 2^/^ reflects the macroscopic mag- 
nification in a Bose-Einstein condensate. To observe the 
MQT, r must be equal to or smaller than the lifetime 
of the BEC, which is of the order of a second. Substi- 
tuting parameters of ^''Rb into Eq. (pT|) . we obtain the 
upper limit for the total number of particles A^max — 36 
for r < 1 s. 

In conclusion, by examining spin-2 BECs we have 
shown that the Bogoliubov theory is, in general, in- 
adequate in identifying the nature and boundaries of 
quantum phase transitions and we should invoke the 
Beliaev theory to resolve this problem. The ground- 
state phase diagram is shown to be significantly af- 
fected by the quantum interference of virtual excitations 
and time-reversal symmetry breaking. This study has 
shed light on the pivotal role of quantum fluctuations 
in the first-order quantum phase transitions. It will 
be an interesting future problem to investigate the im- 
plication of this study to closely related problems such 
as the Coleman- Weinberg mechanism of quantum sym- 
metry breaking [29|, quantum anomaly [soj, and quasi- 
Nambu-Goldstone modes (3l| . 
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